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ABSTRACT 

We report Chandra X-ray constraints for 20 of the 52 high-redshift ultraluminous infrared galaxies (ULIRGs) 
identified in the Spitzer Extragalactic First Look Survey with /„(24//m) > 0.9 mJy, log( ^j 2 ^"*" 1 ) > 1, and 
log(^^^) > 0.5. Notably, decomposition of Spitzer mid-infrared IRS spectra for the entire sample indi- 
cates that they are comprised predominantly of weak-PAH ULIRGs dominated by hot-dust continua, charac- 
teristic of AGN activity. Given their redshifts, they have AGN bolometric luminosities of w 10 45 — 10 47 erg s" 1 
comparable to powerful Quasi-Stellar Objects (QSOs). This, coupled with their high IR-to-optical ratios and 
often significant silicate absorption, strongly argues in favor of these mid-IR objects being heavily obscured 
QSOs. Here we use Chandra observations to further constrain their obscuration. At X-ray energies, we 
marginally detect two ULIRGs, while the rest have only upper limits. Using the IRS-derived 5.8 fim AGN 
continuum luminosity as a proxy for the expected X-ray luminosities, we find that all of the observed sources 
must individually be highly obscured, while X-ray stacking limits on the undetected sources suggest that the 
majority, if not all, are likely to be at least mildly Compton-thick (/Vh>10 24 cm" 2 ). With a space density of 
w(1.4)x 10~ 7 Mpc~ 3 at z^2, such objects imply an obscured AGN fraction (i.e., the ratio of AGN above and be- 
low N H = 1 22 cm" 2 ) of > 1 .7 : 1 even among luminous QSOs. Given that we do not correct for mid-IR extinction 
effects and that our ULIRG selection is by no means complete for obscured AGN, we regard our constraints 
as a lower limit to the true obscured fraction among QSOs at this epoch. Our findings, which are based on 
extensive multi-wavelength constraints including Spitzer IRS spectra, should aid in the interpretation of simi- 
lar objects from larger or deeper mid-IR surveys, where considerable uncertainty about the source properties 
remains and comparable follow-up is not yet feasible. 

Subject headings: galaxies: active galaxies: high-redshift infrared: galaxies X-rays: galaxies 
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1. INTRODUCTION 

One of the most fundamental discoveries borne out of the 
past two decades is that 1) local, massive galaxies almost 
universally host central, super-massive black holes (SMBH; 
> 10 6 Mq) and 2) there exists a tight, fundamental corre- 
lation between a galaxy's SM BH mass, stellar velocity dis- 
persion, and bulge luminosity dKormendv & Richstonel l 19951: 
iMagorrian et alj|1998t iFerrarese & Merrittl 120001) . This key 
result implies that the evolution of a galaxy and its SMBH 
are strongly coupled from a relatively early age. Indirect ver- 
ification comes from the broad agreement between the local 
SMBH mass density, as inferred from the luminosity func- 
tion of bulges (which itself should be intimately related to 
galaxy evolution), and the integrated Active Galactic Nu- 
clei (AGN) emissivity estimated from measurements of the 
Cosmic X-ray Ba ckground (CXRB) and the X-ray luminos- 
ity function (e . g., lYu & T remaine 2002; Marconi et al ]|2004t 
iMerlonil 120041 lLa Franca et alj 120051; Ishankar et al.1 120071) . 
However, direct observations of such symbiotic growth have 
been difficult to obtain. A vital step is to obtain an accu- 
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rate census of all AGN activity at all epochs to truly trace 
the growth of SMBHs. 

Locally, highly obscured (i.e., Compton-thick with 
Nh>10 24 cm" 2 ) a ccretion appears to occur in ~50% of iden- 
tified AGN (e.g.. iRisaliti et all 1 19991 iGuainazzi et"aTll2005t 
Maliz ia et al.l l2009). Similarly strong constraints at high red- 
shift, however, remain elusive. Notably, deep X-ray sur- 
veys have revealed an AGN sky density of > 7200 deg~ 2 , 
at least 2-10 times more than are found at other wave- 
lengths (e.g., iBauer et all 120041) . The bulk of these faint 
AGN have redshifts of 0.5-2.0 and do inde ed appear to 
be ob scured by sign ificant gas and dust (e.g., iBarger et al.1 
l200l ISzokoly etal1l2004t see iBrandt & Hasmgedl2005f for 
a review). Somewhat surprising, however, is the fact 
that such sensi tive X-ray surveys (lAlexander et al.1 [2003; 
iLuo et all 120081) have yet to uncover more than a hand- 
ful of potential Com pton-thick sources similar to the one s 
found locally (e.g., Bauer et all 120041; iTozzi et al l 120061) . 
despite several tentative lines of evidence which suggest 
that Compton-thic k AGN could be at least as plentiful as 
found locally (e.g..lWorslev et al.l20 05; Treis ter & Urry|2005T; 
l Alonso-Herrero et al.l 120061 ; iDonlev et alj 120071: iGilli et al.l 
120071: iFiore et al.ll2.009b . Thus constraints on this hidden pop- 
ulation of high-z Compton-thick AGN must rely on other 
means. 

The mid-infrared (mid-IR) regime offers much poten- 
tial for discovery, since any primary AGN continuum 
that is absorbed must ultimately come out at these wave- 
lengths. Indeed, ultraluminous infrared galaxies (ULIRGs, 
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TABLE 1 

IR-Bright Candidate Compton-ThickULIRGs 



Object 


z 


L 5.8 M m L 3-1000 M m 


J° 

^2-lOkcV 


roc 

L 2-10keV 


Nh 


He 

^2-lOkeV 


References 


Mrk231 


0.0422 


45.10 


46.14 


42.46 


42.46 


24.3 


43.65^44.30* 


1 


Mrk 273 


0.0378 


43.80 


45.78 


42.26 


42.85 


24.2 


43.40* 


2,3 


Mrk 463 


0.0504 


44.80 


45.30 


42.55 


43.36 


23.9 


44.80 


4 


UGC5101 


0.0394 


44.10 


45.59 


42.34 


42.84 


24.1 


44.30 


5,6 


NGC 6240 


0.0245 


43.50 


45.44 


42.09 


42.54 


24.3 


44.20* 


7, 8, 9 


IRAS 00182-7112 


0.3270 


45.70 


46.98 


43.72 


43.90 


24.6 


45.00 


10 


IRAS 09104+4109 


0.4420 


46.10 


46.82 


44.15 


44.15 


24.5 


46.10* 


11 


IRAS 12514+1027 


0.3000 


45.60 


46.57 


42.73 


42.73 


>24.0 


44.20 


14 


IRAS F15307+3252 


0.9257 


46.00 


47.26 


43.78 


43.78 


5^24.0 


45.30 


12 


IRAS 19254-7245 


0.0617 


44.50 


45.70 


42.40 


42.40 


24.5 


44.40* 


13 



NOTE. — Col. (1) Object. Col. (2) Redshift. Col. (3) Logarithm of the rest-frame 5.8 /^m continuum luminosity, calculated 
following S07, in units of erg s~ . Col. (4) Logarithm of the rest-frame 3-1000,/^m continuum luminosity (Lir), in units of erg s~' . Col. 
(5) Logarithm of the observed rest-frame 2-10 keV luminosity of the AGN component derived directly from 2-10 keV data alone, with 
no correction for absorption, in units of erg s~ . Since the X-ray spectra of every source is contaminated by some degree of vigorous 
circumnuclear star formation, the values quoted here should be considered approximate. Col. (6) Logarithm of the rest-frame 2-10 keV 
luminosity of the AGN derived directly from spectral fits to 2-10 keV data alone, corrected for obvious absorption detected only in 
the 2-10 keV band itself, in units of erg s _1 . In some cases, contamination from vigorous circumnuclear star formation precluded any 
assessment of apparent obscuration, in which case the values here are identical to those in Col. (5). Col. (7) Logarithm of the absorption 
column density toward the AGN derived from spectral fitting to all available X-ray data, in units of cm ". Col (8) Logarithm of the rest- 
frame intrinsic 2-10 keV luminosity of the AGN derived from spectral fits to all available X-ray data, assuming the 2-10 keV data is due 
entirely to scattering and/or reflection components, in units of erg s~' . Sources denoted by ' *' have been detected in the 10^4-0 keV band, 
assumed to be largely direct continuum, and should therefore be considered more robust. Nonetheless, all absorption corrections are still 
model dependent and should be considered approximate. Details regarding the X-ray flu x, measured colum n de nsity, and absorption 
corr ections can be fou nd in the references provided in C ol. (9). Col. (9) R eferen ces. 1 — [B raito et al. 2004; 2 — Balestra et al.| 2005t 
3 — Tens et al. 2009; 4 — Bianchi et al. 2008; 5 — Imanishi et al. 2003; 6 — Delia Ceca et al. 2008; 7 — Vianati et al. 1999 8 — 
IKomossa etaU2003t 9 — llwasawa et al. 2009; 10 — Nandra & Iwasawa 20Q3 1 1 — llwasawa etaU200 it 12 — IWilman et aU2003t 
13 — llwasawa et alJ2005t 14 - IBraito et aU2009i 



£ir^10 12 L©) and their low-luminosity brethren have long 
stood out as c andidates for highly obscured accretion onto 
SMBHs (e.g.. ISanders & Mkabej Il996t iFarrah et al.ll2003l 
120071: Nardini et al. 2008), and several of the closest and/or 
brightest members of this c lass have now been confirmed as 
Comp ton-thick AGN (e.g., IComastrill2004t iDella Ceca et alJ 
I2008t see also Table Q). As a class, ULIRGs are al- 
most universally X-ray faint due to their obsc u red nature 
(e^ lFranceschini et alJ 120031: iTeng et all 120051: [Ruiz et alJ 
2007), and have relative AGN contributions which scale 
with increasing infrared (3-1000/im) luminosity, often dom- 
inating the bolometric output in the most luminous ob- 
jects. This population has comparab le space densities 
to Quasi-stellar Objects (QSOs; e.g., ISmail et alJ 119971: 
iGenzel & Cesar skv 2000) and has been proposed a s an impor- 
tant e arly evolutionary stage of these AGN (e.g., iPage etaT] 
12001: lAlexanderetafll2005bt IStevens e t al. 2005 1 ). As such, 
ULIRGs could conceivably host a sizable fraction of the 
most highly obscured, powerful AGN (L DO i ^ 10 46 erg s" 1 ) 
that have evaded detection even in the deepest X-ray surveys. 
The combination of mid-IR color selection, which traces the 
hot dust that obscures AGN at mo st other wavelengths (e.g., 
lLacv et al.ll2004tlStern et al.ll2005l) . and sensitive Spitzer data, 
for instance, routinely offer efficient selection of bright ob- 
scured and unobscured AGN by the hundreds now. Like- 
wise, simple 24/im flux selection may preferentially single 
out AGN-dominated sources in the mid-IR, even when the 
bolometric emis sion from such objects may still be sta rburst- 
dominated (e.g jBrand et alJl2006l:IWatabe et al.ll2009h . 

We focus here on the characterization of a sample of high- 
redshift, 5^/fzer-selected ULIRGs which appear to host very 
powerful obscured AGN based on their mid-IR spectra and 
UV-to-radio spectral energy distributions (SEDs). Our goal 
is to understand how these sources, which represent the ex- 
treme of the overall IR-emitting galaxy population, fit into 
the picture of AGN demography. This paper is organized 
as follows: |j2] describes our high-redshift sample; S}3] details 



our data and reduction methods; fj4] compares our new X-ray 
constraints to existing properties of IR-Bright ULIRGs, ex- 
amines oth er AGN and star forma tion indicators from pro- 
vided by ISaiina et alJ (120071 120081 hereafter S07 and S08, 
respectfully), and investigates how these sources relate to 
other selection techniques; and finally |J5] summarizes our 
findings. We adopt a flat S1a,J1m = 0.7,0.3 cosmology with 
h = // (km s" 1 Mpc _1 )/100 = 0.70, and a neutral hydro- 
gen column density of A^h = (2.5 ± 0.3) x 10 20 cm" 2 for the 
«3.7 d eg~ 2 Spitzer Extragalactic F irst Look Survey (xFLS) 7 
region (Lockma n & Condonl20 05). 

2. SAMPLE 

Our parent sample is comprised of 52 sources 
with f v (2Apm) > 0.9 mJy, log( ^ j 2 "^" ) > 1, and 

log( "^ffl > 0.5 (IYanetal.1 120051) . s elected fr o m the 
full 24fim xFLS catalog. As noted in lYan et all d2007l 
hereafter Y07), this technique selects 59 objects over the 
3.7 deg" 1 xFLS region from the final mid-IR catalogs, of 
which 52 were originally followed-up with Spitzer IRS 
using the initial mid-IR catalogs. Although some objects 
technically lie just below the ULIRG luminosity cutoff, 
for simplicity we hereafter refer to this sample as xFLS 
ULIRGs. The two color criteria are designed to target z ~ 2 
IR -luminous obscur ed objects. As demonstrated in Fig. 2 
of lYan et alJ d2004l) . typical starburst and obscured AGN 
spectral templates migrate into this general color-color region 
by z ~ 1 .5-2, while unobscured AGNs, by contrast, remain 
consistently blue as a function of redshift, and hence are not 
selected. The adopted 24 /im flux density cutoff additionally 
limits our sample to only the most luminous ULIRGs at 
z~2. 

The effectiveness of the color selection is confirmed by low 
resolution, mid-IR spectra (14^1-O^m) taken with Spitzer In- 

7 |http : //ssc . spitzer . caltech . edu/f ls/| 



X-ray Constraints on AGN Properties in z ~ 2 ULIRGS 
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FIG. 1. — Rest-frame Spitzer IRS spectra (solid black curves) and ±1<t errors (grey regions) in units of 1CT 15 W m~ 2 , adapted from Fig. 3 of Sai ina et al.l 
12007). The red diamond shows the 5.8/^m continuum flux extracted from spectral deconvolutio n, which for z < 1 .5 objects had to be extrapolated based on 
additional IRAC photometry (not shown). Also shown is the local starburst galaxy template from Brandl et al. (2006), normalized to the rest-frame IRS data of 
the xFLS ULIRGs to demonstrate the maximum contribution a typical starburst could likely make to the spectra. Note that the templates shown do not constitute 
the actual starburst contributions derived from mid-IR or full SED deconvolutions, which are estimated to be substantially lower in nearly all cases t Sa iina et al.l 
120071 [2"0 08, and further work in preparation). Importantly, while the strong-PAH sources show obvious PAH features comparable to the starburst template, the 
weak-PAH and featureless power-law sources are clearly continuum dominated and look nothing like PAH-dominated starbursts. In nearly all the latter ULIRGs, 
there is strong, excess continuum at 5.8/mi presumably powered by obscured AGN. As noted in blue for individual objects, the case for AGN dominance is not 
as strong in the far-IR, where nearly half (9/20) of the sample are starburst-dominated. 



fraRed Spectrograph (IRS) and optical/near-IR spectra taken 
with Keck. A detailed description of the mid-IR spec- 
tra and analyses combining spectra with far-inf rared and 
sub-m illimeter photometry have been published in I Yan et"aT] 
(120071 hereafter Y07), S07, and S08. Here we summarize the 
salient results for this sample. 

The mid-IR low-resolution spectra provide redshift mea- 
surements for 47 of the 52 sources: the majority 
(35/47 = 74%) lie at 1.5<z<3.2, while a small fraction 
(12/47 = 26%) lie at 0.65 <z< 1.5. One additional source, 
MIPS 279, has a Keck redshift of z = 0.95. The remaining 
four sources are all optically faint (R > 25) and have smooth 
power-law mid-IR spectra: near-IR photometry allows us to 



estimate z p i, = 1.9 ±0.5 for MIPS 15929 based on a spectral 
energy distribution (SED) template fit, while the others are 
only detected in 1-2 bands and thus likely l ie at z p h > 1.5 
(e.g. jAlexander et alJl200ltlRigbv et al.ll2005l) . 

At such redshifts, most of the mid-IR spectra cover the rest- 
frame ~ 5 - 15 /.im band, with a minimum signal-to-noise per 
pixel of 4; for convenience , Fig. Q] reproduces these spec- 
tra from iSaiina et alj d2007l) for all of the xFLS ULIRGs 
listed in Table [2] Locally, the PAH strength has been shown 
to be a robust tracer of star formation (e.g., iBrandl et alj 
l2006t lO'Dowd et alj|2009l) and low-equivalent width sources 
are virtually all AGN-dominated (at lea s t in the mid-IR ; 
e.g., IGenzel et alJ [1991 ILutz et alJ \l99& ITran et al.l l200ll: 
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TABLE 2 
XFLS ULIRGs 



ID 


Net Counts 


Observed Fx 




z 


Rest-frame Flux 


Rest-frame Luminosity 


PAH 


Optical 


T"9.7/jm 




SB 


HB 


FB 


SB 


HB 


FB 




2-10 keV 


2-10 keV 


5.8 fj,m 


3-1000 fim 




AGN 


Marginal X-ray Detections 


8268 


0.6 


5.8* 


6.4 ,v 


<1.60 


6.82 


3.76 


0.80 


<7.24 


<43.34 


44.40 


45.0 


w 




0.8±0.4 


8342 


0.9 


3.6 


4.5"' 


<1.37 


<4.77 


1.75 


1.56 


<1.71 


<43.44 


45.08 


46.0 


w 


yes 


0.2±0.4 


X-ray Upper Limits 


78 


1.2 


3.6 


4.8 


< 1.76 


<3.58 


<5.25 


2.65 


<2.20 


<44.10 


46.24 


46.7 


w 




3.4±0.5 


1 10 


-0.4 


3.5 


3.2 


<1.06 


<5.29 


<5.34 


1.05 


<2.58 


<43.18 


44.79 


45.6 


w 


yes 


0.5±0.5 


1 DJ 


— 1 .z 


— 4.Z 


? ? 

—j.j 


< J. tz 


<lo.ju 


aq 


n on 
u.yu 




^ A1 JA 


AA ^ 
44. JJ 


1 

4 J. J 


w 


yes 


n 04- n c 


279 


-0.1 


0.8 


0.8 


<0.89 


<4.40 


<2.43 


0.95* 


<1.59 


<42.90 


44.50 


45.2 


w 


no 


0.0±0.1 


283 


1.7 


4.1 


5.8 


<1.89 


<9.49 


<4.57 


0.94 


<3.13 


<43.15 


44.43 


44.9 


s 


no 


1.5±0.9 


506 


0.1 


0.2 


0.3 


<1.50 


<5.70 


<3.14 


2.47 


<1.88 


<43.96 


45.40 


46.3 


s 


no 


>6.7 


8135 


1.0 


-0.1 


0.9 


<2.10 


<4.90 


<2.70 


>1.50 p 


<3.76 


<43.73 


>45.25 


>46.0 


w 






8196 


0.0 


-0.2 


-0.2 


<0.86 


<2.85 


<1.57 


2.59 


<1.08 


<43.77 


46.04 


46.5 


w 


yes 


1.3±0.4 


8207 


-1.7 


0.2 


-1.5 


<1.13 


<13.13 


<4.12 


0.83 


<2.44 


<42.92 


44.29 


45.0 


s 


no 


0.9±0.5 


8242 


1.5 


1.5 


3.0 


<1.86 


<7.14 


<3.93 


2.45 


<3.15 


<44.17 


45.65 


46.3 


w 




0.9±0.5 


8245 


0.0 


-0.1 


-0.2 


<0.88 


<5.85 


<1.61 


2.70 


<1.11 


<43.82 


45.81 


46.4 


w 




2.8±0.7 


8493 


-0.7 


-0.4 


-1.1 


<1.16 


<4.80 


<3.17 


1.80 


<3.24 


<43.86 


44.96 


45.4 


s 




>4.7 


15840 


2.6 


0.2 


2.9 


<2.14 


<3.42 


<2.35 


2.30 


<3.18 


<44.11 


45.86 


46.3 


w 




0.2±0.2 


15929 


0.0 


0.9 


0.8 


<0.97 


<4.80 


<2.64 


>1.50* 


<2.46 


<43.54 


>45.22 


>45.9 


w 






15958 


0.9 


0.9 


1.7 


<1.35 


<4.70 


<2.59 


1.97 


<3.07 


<43.93 


45.53 


46.1 


w 




0.5±0.4 


15977 


3.9 


-3.4 


0.5 


<2.36 


<2.53 


<3.48 


1.85 


<3.33 


<43.90 


45.40 


46.2 


w 




0.1 ±0.1 


16006 


-0.8 


0.9 


0.1 


<1.14 


<9.43 


<4.16 


1.90 p 


<2.04 


<43.72 


45.16 


45.9 


w 






16080 


-0.2 


-0.1 


-0.3 


<0.89 


<5.91 


<1.63 


2.01 


<1.60 


<43.69 


45.51 


46.2 


w 


yest 


2.1±0.5 
















X-ray 


Stacking Results 














0.8<z<1.5 


-2.1 


4.5 


2.4 


<0.26 


<3.15 


<1.53 


(0.96) 


<0.65 


<42.44 


(44.53) 


(45.3) 








1.5<z<3.0 


4.3 


3.6 


8.0 


<0.22 


<0.67 


<0.46 


(2.09) 


<0.28 


<42.95 


(45.49) 


(46.1) 









NOTE. — Col. (1) MIPS source number, from Yan et al. 1 2007). Last two entries denote stacking results (see ^3}. Cols. (2^4) Total background-subtracted counts in the 0.5-2 keV 
(SB), 2-8 keV (HB), and 0.5-8.0 keV (FB) bands, respectively, as measured by AOS EXTRACT in 90% encircled energy region at 1.49 keV. 'w' denotes a source detection with 
WAVDETECT at the 10~ 5 significance threshold. Cols. (5-7) Flux in the 0.5-2 keV, 2-8 keV, and 0.5-8.0 keV bands, respective ly, in units of 10~ 15 erg s _I cm -2 assuming a power-law 
model with r = 1.4, If the source is undetected, we instead provide the 3tr upper limit using the Bayesian determination of Kraft et al. i 1991). Col. (8) Redshift, from S08. "p" 
indicates a photometric redshift. Note that MIPS 279 ("*") was quoted in early xFLS papers as having a redshift of 1 .23; this as since been revised to the value adopted in the table. Col 
(9) Observed rest-frame 2-10 keV flux in units of 10~ 15 erg s -1 cm -2 assuming a power-law model with r = 1.4. Col (10) Logarithm of the observed rest-frame 2-10 keV luminosity 
in units of erg s -1 . Col. (11) Logarithm of the rest-frame 5.8 fim continuum luminosity, from S07, in units of erg s -1 . Co! (12) Logarithm of the rest-frame 3-1000 fim continuum 
luminosity, from S08, in units of erg s -1 . Col (13) Relative strength of PAH features in IRS spectrum, from S07. Strong ("s") PAH sources have EW(7.7/im)>0.8/i.m, while weak 
("w") PAH sources have lower EWs. Col. (14) Presence of AGN as assessed by optical spect roscopy, from S08. M IPS 16080 ("^") has starburst-like line ratios, but an asymmetric 
[O III] profile typical of an AGN outflow; this is similar to brighter, radio-excess ULIRGs (e.g., Buchanan et al. 2006). Col. (15) Silicate absorption at 9.1 fim in dimensionless units. 



iFarrah et alj 120071: iDesai et alj 120071) . Despite the modest 
signal-to-noise for some spectra, comparison to m aximal con- 
tributions from an average starburst template of Brandl et al. 
(2006) in Fig. Q] clearly demonstrates that the vast majority 
of xFLS ULIRGs studied here are continuum-dominated. We 
highlight that aside from the four starburst-dominated objects, 
the best-fitted starburst contributions to the mid-IR are gener- 
ally much lower than the maximal contribution shown. 

When the entire xFLS sample are more rigorously decom- 
posed into PAH, continuum, and obscuration spectral com- 
ponents, S07 find that ^75% of the sample are weak-PAH 
sources with £W / - fs ,(7.7/.tmPAH) < 0.8 /im, consistent with 
AGN-dominance in the mid-IR. Although the remaining 25% 
appear to be starburst-dominated, the S07 template decon- 
volution indicates that only ^50% of the 5.8/im continuum 
can reasonably be attributed to star formation, prompting 
the need for an additional hot-dust component from a dust- 
obscured AGN. These results are reinforced by the known 
correlation between star burst PAH strength and the host 
galaxy stellar bumps ( e .g.jLacv et a l. 2004; Stern et aLll2005b 
IWeedman et alj 120061; iTeplitz et alj 12007b . wherein strong- 
PAH sources generally show clear stellar bumps in the IRAC 
bands while weak-PAH sources show no such bumps and usu- 
ally are just power laws (S07). Subsequent analysis of the op- 
tical spectra, radio properties, and UV-to-radio SEDs by S08 
further confirm the mid-IR spectral deconvolutions, conclud- 
ing that strong-PAH sources typically have AGN contribu- 
tions of ^20-30% of the total Lir, while weak-PAH sources 
generally have AGN contributions of >10%. Notably, weak- 



PAH sources are roughly twice as likely to lie at z > 1.5 
and have a substantially higher fraction of optically-identified 
AGN compared to strong-PAH sources (see Fig. Q] and S08). 
Thus the conclusions regarding z~2 ULIRGs should be ro- 
bust. 

We additionally note that Sajina et al. (2010, in prep.) study 
the UV-to-radio SEDs of a large sample of 24um-selected 
sources including most of the objects discussed here. This 
work includes comparison with local galaxy templates as well 
as fits to composite quasar+starburst templates. They find 
that such fits are not sensitive to the details of the IRS spec- 
tra but rather to the overall 1-1000/im SED shapes of these 
sources. They conclude that although weak-PAH sources 
can still have their total IR luminosities dominated by star- 
burst activity, the 5.8/im continuum of such sources is over- 
whelmingly dominated by obsc ured quasars (e.g., Fig. [T]). 
This conclusion supports those o flPolletta et all ((2008), whose 
sample also includes a few of the sources studied here, but 
where the fitting includes theoretical torus models rather than 
a quasar template. Both studies suggest that unless z^2 star- 
bursts and AGN have dramatically different SED shapes than 
has been hitherto observed, the conclusion that weak-PAH 
sources have 5.8/.tm continua dominated by AGN is sound. 

The total IR luminosities of the sample, derived from the 
mid-IR spectrum plus 70/im, 160/im, and 1.2mm photom- 
etry (S08), range from (0.04-2) x 10 13 Lq, with a median 
of 5xlO 12 L0. While of comparable bolometric luminosi- 
ties to QSOs, these sources are typically much more ob- 
scured, as evidenced by their steeper mid-IR slopes and of- 
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ten significant 9.7 /im silicate absorption (e.g., mean optical 
depth (rg.-j) ~ 1.4). This latter quantity, T9 7, appears t o cor- 
relate strongly with X-ray derived N H dShi et al.l 120061) . We 
note that the distribution of slope and absorption values are 
some what more extreme than those ty pical of type II AGN 
(e.g JDesai et al.l2007tlHao et al.l2007l) . suggesting that these 
objects may not merely be QSOs observed at unfavorable 
viewi ng angles (e.g., the AGN unification model; lAntonuccil 
1993), but might represent a particular obscured class or evo- 
lutionary phase of QSOs. 

Clearly there could be variations in the individual starburst 
and AGN mid-IR SEDs, and thus the accuracy of the spec- 
tral decompositions is somewhat limited. However, this is un- 
likely to change the majority of 5.8/im continuum measure- 
ments by more than 10-20% typically, and thus would have 
little effect on our results. We also caution that the quoted 
rest-frame mid-IR luminosities from S07 have not yet been 
corrected for intrinsic mid-IR absorption, which is difficult to 
determine on a source-by-source basis. In particular, some 
fraction of our sources have red continua indicative of strong 
mid-IR extinction but no obvious accompanying silicate ab- 
sorption, and thus may have considerable variation in either 
their dust geometry or intrinsic extinction law. Although our 
quoted continuum at 5.8 urn lies near a minimum in the mid- 
IR extinction curve (e.g., Chiar & Tielens 2006), the impact 
of dust could still be considerable. With extinction typical 
of the Galactic Center (GC), for instance, our T9.7 measure- 
ments could equate to a 5.8 fim flux decrement of — 1.6 on 
average and — 10 in a few extreme cases. Additionally, the 
presence of molecular absorption features due to water ice 
and hydrocarbons often accompany strong silicate absorp- 
tion and may result in additional intrinsic abso rption (e.g., 
around 6.2/zm) in some weak-PAH sources (e.g..lSpoon et al.l 
2001 12001 llmanishi et alJl2(M l2(M iRisahti et alJlSOQft 
Sani et al. 2008). Care has been taken to exclude affected 
spectral regions but such features m ay still impact the decom- 
position fitting. ISaiina et"afl (120091) explores in detail some 
instances of water ice and hydrocarbons within the xFLS sam- 
ple, to which we refer interested readers. Both features could 
lead to significant underestimatates of the AGN continuum, 
which must thus be regarded as lower limits to their intrinsic 
values. Finally, objects with z < 1.5 have somewhat larger 
uncertainties in Ls.s^m because estimates can only be extrap- 
olated based on IRS spectral decomposition above rest-frame 
— 6-8/im and broadband IRAC photometry below. 

Thus to summarize the strong evidence for dominant AGN 
activity in the xFLS sample: 

• Based on spectral deconvolution of the IRS spectra, 
k.15% of the sample are continuum-dominated, weak- 
PAH objects powered by AGN in the mid-IR, and even 
among strong-PAH objects, additional AGN continuum 
components are generally required. While there clearly 
could be some uncertainty in the deconvolution, the 
possible issues (e.g., templates, extinction from dust 
and ice) largely skew toward underestimating the true 

continuum of the xFLS ULIRGs. 

• Spectral decomposition of the rest-frame UV-to-radio 
SEDs confirm the modeling results of the mid-IR data 
alone, suggesting that strong-PAH sources typically 
have -20-30% AGN fractions of L m , while weak-PAH 
sources by contrast tend to have >70% AGN fractions, 
with a few outliers having comparable contributions 



from AGN and starbursts. 

• Radio and, where available, optical-line diagnostics 
support the presence of AGN in —60% of the z>1.5, 
predominantly weak-PAH ULIRGs, independent of any 
IR-based diagnostics. 

To explore the X-ray energetics and constrain the nature 
of any potential X-ray absorption in the above sources, we 
obtained X-ray observations for a random subset of the full 
sample. 

3. OBSERVATIONAL DATA AND REDUCTION METHODS 

We observed 20 of the 52 xFLS ULIRGs with five Chandra 
ACIS-I pointings of 30 ks each (PI: Yan; Obsids: 7824, 7825, 
7826, 7827, and 7828). The data were processed following 
standard procedures using CIAO (v3.4) software. 8 We addi- 
tionally removed the 0"5 pixel randomization, corrected for 
charge transfer inefficiency, performed standard ASCA grade 
selection, excluded bad pixels and columns, and screened for 
intervals of anomalous background using a threshold of 3a 
above or below the mean (this excluded at most 1-2 ks for 
each observation). Further analysis was performed on repro- 
cessed Chandra data, using CIAO, FTOOLS (v6.3) and custom 
software including AOS EXTRACT (v3.172) 9 

Source detections were determined using the WAVDETECT 
source-searching algorithm with a 10~ 5 significance thresh- 
old. While this threshold is typically considered 'lenient' 
(i.e., 10-20 false sources over full Chandra ACIS-I field), the 
low xFLS ULIRG source density means we are only inter- 
ested in detections at a handful of locations over the ACIS- 
I field-of-view, and thus it translates here into a very robust 
search significance. Since our sources are likely to be ob- 
scured, we searched in the standard full (0.5-8.0 kev), soft 
(0.5-2.0 keV), and hard (2-8 keV) bands separately. We 
additionally searched in the observed 1-4 keV band, which 
roughly equates to rest-frame 2-10 keV at z > 1.5, to deter- 
mine the assignment of detections and upper limits at the rest- 
frame 2-10 keV fluxes assessed below. There are typically 
—90-100 X-ray sources detected per observation, which we 
matched to optical counterpar ts from the xFLS optical sur- 
vey catalog dFadda et al.l 120041) . providing —30 highly solid 
matches within 0."5 per observation. The matches allowed 
registration of the X-ray astrometric frame to the optical, 
with typical linear shifts of 0."l-0."2 to the original X-ray 
frame. The resulting la registration residuals are «0."3, en- 
suring spatially robust X-ray identifications and upper limits. 
Among the 20 objects observed, only two were formally de- 
tected in any of the above bands (see Table|2]i. 

Aperture-corrected photometry was performed using 
ACIS EXTRACT with a 90% encircled-energy region derived 
from the Chandra PSF librar y, with 3a upper limits calculated 
following iRraft et all (1199 ll) . Fluxes were derived assuming 
a power-law model with T = 1 .4 consistent with the spectrum 
of the CXRB. We also performed photometry on specific sub- 
bands to allow more accurate constraints on rest-frame fluxes 
at the redshifts of our sources. To this end, we use the ob- 
served 0.5-2.0 keV and 1^4 keV bands to estimate the rest- 
frame 2-10 keV fluxes (or upper limits) and subsequent lumi- 
nosities for objects with z < 1 .5 and z > 1 .5, respectively. An 
additional K-correction of order 15-20% is applied to account 
for individual redshifts. 

8 http://asc.harvard.edu. 

9 http://www.astro.psu.edu/xray/docs/TARA/ae_users_guide.html 
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Regarding the two detected sources, we find that one is 
detected in the full and hard bands, while the other only in 
the full band. Although the photon statistics here are poor 
(both sources are formally c onsistent with zero net counts 
at l-2cr confidence following Gehrels 1986), the photon en- 
ergy distributions lean toward both being heavily obscured. 
An unabsorbed AGN with a T = 1.9 power-law spectrum, 
for instance, would have three times as many 0.5-2.0 keV 
counts as 2-8 keV counts, easily detectable in our Chan- 
dra observations. Given the effective area of Chandra, a 
large neutral hydrogen column of ~10 23 — 10 24 cm -2 would 
be required to arrive at the observed count distributions of 
both sources. Neither source is formally detected in the 
rest-frame 2-10 keV bands described above, however, so 
we quote their 3<r upper limits for this flux. Assuming 
T = 1 .4 again (T = 1.9 would change these values by only 
«10%), the observed full-band luminosities of MIPS 8268 
and MIPS 8342 are 1.0 x 10 43 erg s" 1 and 3.2 x 10 43 erg s" 1 , 
respectively, and extrapolate to rest-frame 2-10 keV lumi- 
nosities of 7.8 x 10 42 erg s" 1 and 2.5 x 10 43 erg s" 1 for the 
two objects. 

To place stronger average constraints on our sample, we 
also performed X-ray stacking analyses on the undetected 
sources, which has been successfully employed on numerous 
sourc e populations in both wide and deep field X-ray surveys 
(e.g.jBrandt et al.ll2001l:lBrand et al.ll2005HDaddi et aT]|2007l: 
Lehm er et al.l2008l) . We divided the 20 ULIRGs into two sub- 
sets based on redshift (0.8<z<1.5 and 1.5<z<3.0), taking 
care to exclude both the two detected sources and the three 
sources which lie on ACIS chip S7, since the extended PSFs 
from the latter include too much background to improve the 
signal-to-noise. For consistent aperture corrections, we again 
used the 90% encircled-energy aperture region measured at 
1.49 keV from ACIS EXTRACT. For rest-frame luminosity es- 
timates, we stacked counts only from the 2-10 keV rest-frame 
bandpasses of the sources (i.e., in the observed 0.5-2.0 keV 
and 1.0^4.0 keV bands for z>1.5 and z<1.5, respectively). 
Aside from the use of individual geometric source apertures 
in stead of circul a r aper tures, our method largely follows that 
of iLehmer et al.l (120081 and references therein). To properly 
account for spatial variations in pixel sensitivity due to chip 
gaps, bad pixels, and vignetting, the total number of back- 
ground counts within our stacked aperture was determined by 
scaling the cumulative background counts found within lo- 
cal annuli by the ratio of the summed exposure times in the 
source and background regions, respectively. No significant 
detection was found in either redsh ift subsamp l e, 10 so we cal- 
culated 3(7 upper limits following Kraf t et al.l (11991b - These 
count limits were converted to 2-10 keV rest-frame fluxes 
using the previously adopted spectral model for photometry 
above. 

4. DISCUSSION 

We would like to use the above constraints to understand the 
nature of the xFLS ULIRGs and their context with respect to 
the rest of the high-z AGN population. Continuum emission 
at hard X-ray and mid-infrared wavelengths, as well as emis- 
sion from the narrow-line region, are widely considered to 
provide the m ost robust constraint s on AGN bolome trie lumi- 
nositi es ("e.g..lBassani etalJ 119991: IXu et al.l|1999b iLutz et al l 
I2004t IHeckman et al.l l2005Ulmanishil 120061: IMelendez et all 

10 The most significant detection was only al.l a in the full-band for the 
1.5<z<3.0 subsample. 
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2009), as each offers a unique measure of the primary AGN 
energy output over large dynamic range (^4-5 dex) before 
contamination from host star formation sets in. The depen- 
dence of AGN properties on orientation and intrinsic obscu- 
ration is a long-standing problem when tr ying to assess the 
true power of AGNs (e.g. jAntonuccil 19931) and none of these 
tracers are without its faults. However, because these trac- 
ers are affected differently by orientation and obscuration, 
we can assess the nature of obscuration present in our xFLS 
ULIRGs by comparing their X-ray upper limits to other less 
obscured tracers such as rest-frame 5.8/im continuum and (in 
a few cases ) [O III1 flux. This is es sentially the same approach 
adopted by Alexande r et al.l d2008, hereafter A08) to constrain 
Compton-thick AGN in the GOODS-N region. First, how- 
ever, we would like to assess the applicability of these rela- 
tionships to ULIRGs using several bright, well-characterized 
objects from the literature. 

4.1. Guidance from IR-Bright ULIRGs 

A correlation between 2-10 keV and 5.8/im luminosity has 
now been well-established fo r both type 1 and absorption- 
corre cted type 2 AGN (e.g., ILutz et al.l 120041: ISturm et all 
2006). While it is not yet clear why this correlation is as tight 
as observed (e.g., the apparent absorption-corrected 2-10 keV 
luminosities of type 2 AGN can often be substantially 
lower than extrapola t ions based on >10 ke V luminosities; 
IHeckman et all 120051: IMelendez et alj [2008), or what func- 
tional form best describes the correlation is (e.g.. ILutz et al.l 
l2004tlFiore e7aTll2008l: lEanzuisi et al.ll2009l) . empirically this 
correlation implies that AGN are relatively robust, immutable, 
and scalable physical systems. In both panels of Fig. [2] we 
show the intrinsic rest-frame 2-10 keV-to-5.8 /im range for 
type 1 AGNs (i.e., only continuum-dominated mid-IR sources 
with little or no PAH emission), where the correlation is given 
asL 5 .8 M m = lO^'^^Lj-ioteT' with both luminosities in units 
of erg s , and the grey shaded region denotes la luminosity 
range (Bauer et al., in preparation); we note tha t this fit is 
consis tent with the original correlation found by ILutz et ail 
(2004) at low lu minosities and lie s int ermediate betw e en the 
extrapolations of ILutz et all (12004 and lLanzuisi et all (120091) 
at high mid-IR luminosities. 

To demonstrate how X-ray absorption and mid-IR extinc- 
tion affect the intrinsic AGN output, we also show in Fig. [2] 
the expected ratios for a typical AGN (taken as the center of 
the unobscured AGN range) assuming (1) the 5.8/im contin- 
uum is intrinsic but the X-ray emission is absorbed by a col- 
umn densities of Na ~ 10 24 cm -2 and A^h ~ 10 25 cm" 2 , 11 and 
(2) the X-ray emission is intrinsic but the 5.8/im continuum 
is extincted by A^^ m ~ 2 mag, 12 as might be the case in many 
ULIRGs. Both decrements and their associated errors con- 

1 The value of Na ~ 10 24 cnT 2 is roughly where obscuration is ex- 
pected to start becoming Compton -thick, and has b een ca lculated using a 
model similar to that presented in Alexander et al. (2005a) and IGilli et alj 
12007), wherein the primary radiation is obscured and we only detect re- 
flection and scattering components that typically comprise f«5% of the in- 
trinsic 2-10 keV emission. By Nh ~ 10 25 cm -2 , the source is expected to 
be 'fully' Compton-thick such that no primary continuum escapes, and has 
been calculated assu ming the best-fitted X-ray spectrum of NGC 1068 from 
Bass ani et al] fT999h . 

12 If we comb ine the Cardelli et al] 119891 optical/near-IR and 
Chiar & Tielens 2006 mid-IR extinction curves, A6 Mm ~ 2 mag is equivalent 
to A\ ~ 36 mag assuming Ry = 5 appropriate for dense starbursts/H II 
regions or Ay ~ 42 mag assuming R\ = 3.1 appropriate for Milky Way 
interstellar material, respectively. 
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FIG. 2. — Rest-frame 2-10 keV luminosity versus 5.8 fim AGN continuum luminosity for several bright ULIRGs from published literature thought to host 
Compton-thick or nearly Compton-thick AGN (left) and the Cftanrfra-observed xFLS ULIRGs (right). All 5.8 fim AGN luminosities were calculated from the 
best-fitting AGN component via spectral deconvolution (e.g., S07) and have not been corrected for mid-IR absorption. In both panels, the grey shaded region 
denotes the Icr scatter in rest-frame 2-10 keV-to-5.8 fim luminosity ratio range found for type 1 AGNs (Bauer et al., in preparation), the black dotted lines 
demonstrate the effects of X-ray absorption by column densities of ~ 10 24 cm -2 (using the model presented in Alexander et al. 2005a) and Nu ~ 10 25 cirT 2 
(assuming the X-ray spectrum of NGC 1068; Bassani et al. 1999), and the grey dashed lines show the effect of mid-IR extinction by A^^ m ~ 2 mag to unobscured 
and Nyi ~ 10 24 cirT 2 sources {Chiar & Tielens 2006). For the bright, individually-labeled ULIRGs (left): red lower limits denote absorbed, rest-frame Z.2-l0keV 
values (i.e, the observed data; note the composite spectra for some objects were too contaminated by star formation at low energies to assess the AGN contribution 
or apparent absorption unambiguously); blue stars denote absorption-corrected, rest-frame Z-2-lOkeV values derived from the best-fitted models to < 10 keV 
spectra, assuming the data are predominantly direct continuum; green stars denote absorption-corrected, rest-frame L2-l0keV va ' ues derived from the best- 
fitted models to both < 10 keV and >10 keV spectra, assuming the >10 keV data is largely direct continuum while the < 10 keV data is only reflected/scattered 
continuum; and green upper limits are absorption-corrected, rest-frame Z.2-l0keV extrapolations derived from the best-fitted models to < 10 keV spectra, assuming 
the data is reflected/scattered continuum only. Note that all luminosities were taken from the literature. For the xFLS ULIRGs (right): X-ray upper limits (3cr) 
are calculated following Kraft et al. ( 1991), assuming T = 1.4; blue points indicate WAVDETECT full-band detections; red points denote sources with photometric 
redshift constraints; and green bars represent the average stacked upper limits for our z ~ 1 and z ~ 2 subsamples. The individual upper limits lie close to the 
Na ~ 10 24 cirT 2 track, demonstrating that most xFLS ULIRGs are highly obscured, while the stacked averages imply that a significant fraction are likely to be at 
least mildly Compton-thick. For the five xFLS ULIRGs with observed [O III] constraints, we have converted the [O III] luminosities (uncorrected for extinction) 
to equivalent absorption-corrected X-ray luminosities using the correlation from Bauer et al. (2009, in preparation); those observed by Chandra are plotted as 
grey upper bars, while the rest are shown as crosses. 



tribute to the observed scatter of this relation, and thus it is 
natural to wonder how ULIRGs, which are some of the most 
dusty and obscured objects in the Universe, fit into the above 
picture. 

To this end, we assembled a sample of IR-bright ULIRGs 
from the literature, all of which have high quality Spitzer 
IRS and 2-10 keV spectra and are widely considered to host 
Compton-thick AGN. This list is regrettably small because 
such objects are typically X-ray faint and difficult to con- 
strain. Importantly, half of the sample have published X-ray 
detections above 10 keV with Suzaku or BeppoSAX, allowing 
the best assessment of the apparent direct X-ray continuum to 
date. Basic properties of this sample are listed in TableQ] 

The observed and intrinsic X-ray luminosities of these IR- 
bright ULIRGs were taken directly from the references in Ta- 
ble Q] the X-ray data reduction and spectral fitting rely on 
the assumptions detailed therein. For Mrk 231, we list both 
of the intrinsic X-ray luminosities provided by Brai to et all 
(2004), depending on whether the X-ray spectrum is domi- 
nated by scattering (low value) or reflection (high value). For 
both Mrk 463 and UGC 5101, it has been argued that the low 
observed Fe Ka equivalent widths (^200-400 eV) demon- 
strate that the direct AGN continuum is being observed. If 
scattering dominates over reflection, however, such equivalent 
width s can still be fully consis tent with Compton-thick AGN 
(e.g., iMurphv & Yaqoobl r2009). as in fact appears to be the 



case for Mrk 231. Thus for these two sources, we compute 
intrinsic X-ray luminosities assuming the observed continua 
are reflection/scatter dominated and compr i se ~5% of their 
intrinsic values (e.g., Alexander et al. 2005a; Gilli et ai1l2007l 
hereafter G07; Mol ina et alj|20 09): these corrections are de- 
noted with dashed green upper limits in the left panel of Fig. [2] 
since they are more tentative and open to interpretation. 

The 5.8/im AGN continuum luminosities for this bright 
ULIRG sample were derived following the procedure used 
in S07. Specifically, IRS low-resolution spectra in the form 
of pipeline-processed, Basic Calibrated Data (BCDs) were 
downloaded from the Spitzer archive. Sky backgrounds were 
generated by subtracting BCD images at two different node 
positions. Flux- and wavelength-calibrated 1-D spectra were 
then extracted from the sky-subtracted 2-D, BCD images us- 
ing the SSC spectral extraction software, SPICE (v2.2). 13 Fi- 
nally, four segment spectra covering 7-35/im were are fit to- 
gether to eliminate bad pixels at the edge of each order. The 
resulting spectra were corrected for Galactic extinction and 
fit with a power-law continuum and PAH template of M 82. 
Given that M82 has one of the strongest observed mid-IR 
dust continua, our constraints on AGN hot-dust components 
should be considered relatively conservative. In cases where 
prominent 6pm ice absorption is present, we masked it out 

13 http://ssc.spitzer.caltech.edu/postbcd/spice.html 
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and fit only the continua on either side. The Ls^m values are 
the monochromatic rest-frame luminosities read off from the 
best-fit models. In all cases except NGC6240, the continua 
are dominated by a power-law component. This is also the 
case for our high-z xFLS sample, which is largely composed 
of weak-PAH sources. As discussed in Ej2j determining the 
starburst contribution to Ls.s^m in the few strong-PAH sources 
where it may be significant is complicated by the highly un- 
certain levels of hot dust emission that can be associated with 
pure starbursts. Using M 82 as a template, for instance, we 
estimate that up to half of L5.8 Mm in NGC 6240 could be due 
to its starburst. A starburst contribution of <30% can also be 
seen in Mrk273 and UGC5101, but is negligible in the rest 
of the IR-bright ULIRG sample. 

The observed 2-10 keV luminosities of the IR-Bright 
Compton-thick ULIRGs in the left panel of Fig. [2] all lie 
around the Nn ~ 10 24 cm" 2 track, confirming our expectation 
that they are highly obscured. Once the effects of known 
X-ray absorption are accounted for (green limits), however, 
these Compton-thick ULIRGs generally lie on or above the 
type 1 AGN correlation. Importantly, the case for large X-ray 
corrections is strongest for the sources with solid >10 keV 
X-ray detections (green stars), which have the least amount 
of ambiguity regarding their intrinsic X-ray luminosities. 

Considerable uncertainties in the mid-IR also remain from 
the poorly constrained gas-to-dust ratios and dust geometries 
in ULIRGs. Theoretical arguments suggest that the mid-IR 
emission can perhaps vary by up to an order of magnitude 
depending on the distribution and composition of the ob- 
scurin g dust (e.g., iPier & KrolikHl992t iNenkova et al.ll2002l 
2008). If the gas-to-dust ratios were substantially lower for 
ULIRG AGN, for instance, then the mid-IR luminosity con- 
straints could overestimate the true bolometric power of the 
AGN. Such a correction might move some ULIRGs into bet- 
ter alignment with the typical type 1 range, but would make 
current discrepancies with others even more extreme. Particu- 
lar dust geometries/orientations, on the other hand, could lead 
to a few magnitudes of extinction even at mid-IR wavelengths 
and cause us to underestimate the true mid-IR luminosities 
from these objects. The A(,^ m ~ 2 mag track demonstrates the 
potential degree of such an effect. 

The left panel of Fig. [2] illustrates that the known correla- 
tion between intrinsic rest-frame 5.8/im and 2-10 keV emis- 
sion in AGN does effectively extend to ULIRGs and that ob- 
served 2-10 keV decrements should provide realistic diagnos- 
tics of obscuration in all objects containing AGN. Settling on 
a sensible criteria for selecting Compton-thick AGN, how- 
ever, is not clearcut. A conservative approach is to adopt the 
A^h ~ 10 24 cm" 2 track as our Compton-thick criteria. While 
this neglects the potential mid-IR extinction that is likely to 
accompany the X-ray absorption, it ensures that the selected 
sample of candidates will have negligible contamination from 
less obscured sources. Under this scheme, only six of the 
above 10 widely -regarded Compton-thick IR-bright ULIRGs 
would be selected, demonstrating that we might miss a con- 
siderable fraction of likely Compton-thick AGN with such 
an approach. On the other hand, adopting the more inclu- 
sive A^h ~ 10 24 cm" 2 +A6 M m ~ 2 mag track as our Compton- 
thick criteria would select all of the IR-bright ULIRGs, but 
may potentially select AGN which are somewhat less heavily- 
obscured as well. Clearly there are merits and drawbacks to 
each criteria, so we will employ both in the next section to 
gauge the nature of the xFLS ULIRGs. 



4.2. Constraints on the Nature ofxFLS ULIRGs 

In the right panel of Fig. [2] we compare rest-frame 5.8 /im 
AGN continuum and 2-10 keV luminosity constraints for 
xFLS ULIRGs. For a given mid-IR luminosity, we typi- 
cally find a large decrement between our individual X-ray up- 
per limits and the expected 2-10 keV emission, as estimated 
from the intrinsic X-ray-to-mid-IR correlation, demonstrat- 
ing that all of the xFLS ULIRGs should be heavily obscured 
AGN. Only two of the individual 3cr upper limits lie be- 
low the conservative A^h ~ 10 24 cm" 2 track, while all but two 
lie below the inclusive A^h ~ 10 24 cm~ 2 +A(,^ m ~ 2 mag track. 
Thus some fraction of xFLS ULIRGs are likely to be at least 
mildly Compton-thick, but many individual upper limits leave 
ambiguous constraints. 

The locations of the two marginally full-band detected 
xFLS ULIRGs, MIPS 8268 and MIPS 8342, can provide some 
initial guidance. As noted in Sj3) neither are formally detected 
in the rest-frame 2-10 keV sub-band, and thus appear in Fig. [2] 
as upper limits. The rest-frame 2-10 keV luminosities extrap- 
olated from their observed full-band counts lie below their 
upper limits by factors of 2.8 and 1.1, respectively, but still 
roughly a factor of 2 above the conservative Nn ~ 10 24 cm" 2 
track. Both extrapolations are at least broadly consistent with 
our crude band ratio analysis in ij3] If these two objects repre- 
sent the least obscured xFLS ULIRGs, then we would expect 
the rest to lie near or below the conservative track accordingly. 

We can improve upon these constraints by stacking the 
X-ray counts of the undetected sources in two subsamples, 
split into crudely matching redshift ranges (see Sj3j. Impor- 
tantly, we find that neither sample is statistically detected, 
with stacked 3a upper limits lying factors of ml. 6 and «5.2 
below the conservative A^h ~ 10 24 cm" 2 track for the z ~ 1 
and z ~ 2 samples, respectively. Assuming all sources con- 
tribute roughly equally to the stacked signal, these low lim- 
its imply that >70% of the z ~ 1 and >90% of the z ~ 2 
C/ianc/ra-observed xFLS ULIRGs should be at least mildly 
Compton-thick at 99% confidence, and perhaps some frac- 
tion even fully Compton-thick (A^h > 10 25 cm" 2 ). Assuming 
the C/iant/ra-observed objects are representative and follow a 
normal distribution, then we can calculate 'margin of error' 
constraints 14 on the potential Compton-thick fraction for the 
entire sample. For z ~ 1, we only stack 3 of 15 total xFLS 
ULIRGs, resulting in a candidate Compton-thick fraction of 
>25% at 90% confidence, but no meaningful constraint at 
higher confidences. For z ~ 2, we stack 11 of 37 total xFLS 
ULIRGs, resulting in candidate Compton-thick fractions of 
>80% and >70% for confidences of 90% and 99%, respec- 
tively. Our general conclusion that a majority, if not all, xFLS 
ULIRGs are likely to host luminous Compton-thick AGN is 
supported by the fact that the individual xFLS upper limits lie 
in the same range as some observed X-ray lu mino sities from 
our IR-bright Compton-thick ULIRG sample (i H.ll l. implying 
that the xFLS objects could be even more extreme. 

This result is predicated on the fact that the X-ray decre- 

14 When only a subset of the total population is sampled, the associated 
error (often termed 'margin of error') on the fraction of the total population 

with a particular property can be calculated using z\J p ^~ p ^ "/, where z is the 

'critical value' related to the form of the underlying distribution, the degrees 
of freedom, and the confidence interval of interest (in this case a 1 -sided t- 
distribution), p is the fraction of a particular property found from the subset 

population, and / is the finite population correction fact or W , n is th e 
subset number, and N is the total population number (e.g., Peck et al. 2008). 
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ment provides a reliable estimate of the obscuration, when 
compared to t he AG N continuum luminosity at 5.8 /im. While 
we argued in 8 14.1 I that this procedure is robust, emission-line 
luminosities from AGN can provide an independent check on 
the AGN power estimated from our 5.8 measurements. 
Only five xFLS ULIRGs have robust [O III] detections among 
nine observed (S08), which we plot in the right panel of Fig. [2] 
These emission-line luminosities have not been corrected for 
contamination by star formation or extinction, primarily be- 
cause neither is well quantified. The expected contribution 
to the observed L[oni] from star formation should be minimal 
(< 1—10%), as long as the well-established star- formation cor- 
relations between [QUI], [O il], andL/« (e.g.. Ho pkins et alj 
l200llMoustakas et alj|2006h hold for the xFLS ULIRGs. As 
discussed in S08, the effects of internal extinction are more 
problematic, since the few sources with direct rest-frame op- 
tical constraints exhibit a wide range of apparent extinctions. 

For instance, among the five [O III] -detected xFLS 
ULIRGs in Fig. |2j direct assessment of the optical extinc- 
tion via the Balmer decrement is available for four. Assum- 
ing no systematic errors exist between near-IR spectral seg- 
ments, two sources show no evidence for extinction, while 
the other two have E(B-V) = 1.2 and E(B-V) > 2.2, re- 
spectively. Among comparabl e IR-bright, optically-faint z ~ 
2 sources, Brand et al. (2008) also report extinctions in the 
range E(B-V) = 1.0-1.9. Thus the amount of optical extinc- 
tion can vary dramatically from source to source and poten- 
tially be quite severe. There also remains some possibility that 
the [O III] emission does not see the same extinction as the 
Balmer lines. We can alternatively estimate E(B-V) based on 
the silicate absorption in the mid-IR, which lies in the range 
T9.7 ~ 1 .3-2.7 for the five sources and provides constraints on 
extinction from the nuclear region itself. The ratio Ay / 79.7 
is principally dependent on the relative abundances of sili- 
cate and graphite grains: «18 for diffuse interstellar material 
(ISM), ^ 18-40 for denser clouds, and &9 near the GC (e.g., 
iRoche & AitkeiJl985tlWhIttell2^0ltlChTar et al.ll2007h . If we 
conservatively adopt the low GC conversion and Ry = 5, then 
we would expect Ay ~ 12-24 or E(B-V) ~ 2.3-4.9, which 
is substantially higher than the direct estimates. Even adopt- 
ing a somewhat conservative E(B-V) ~ 1 .0, the correction to 
[OIII] is already ~l-2 dex. 

The five [O III] -detected xFLS ULIRGs have 
log[L[oiii]/V£i/(5.8^iiri)] ratios of -3.0 to -3.5, which 
systematically lie 1-2 dex lower than the a verage type 1 and 
type 2 AGN/ULIRGs, respectively (e.g., lHaas et all 120071; 
Bauer et al., in preparation). The discrepancy is consistent 
with our expectation that the [O III] luminosities are highly 
extincted. Thus our uncorrected values should be considered 
very conservative lower limits to the true AGN power. 

To compare against our X-ray limits in the right panel 
of Fig. [2] we convert the [OIII] constraints to X-ray 
ones using the [O III] -to-X-ray correlation for type 1 AGN 

(L[onrj = 10 6 ' 9±2 ' 2 L2_ioteV° 5 )' where the emission-line lumi- 
nosities have not been corrected for extinction (Bauer et al., 
in preparation); both luminosities are in units of erg s" 1 . 15 

15 We note that this £ [oiin/£2-l0ke V corre l ation differs significantly 
from those found by , e.g.. INetzer et alj 0.006), Panessa et al. (2006) and 
Melendez et al. ( 2008), which primarily appears to be due to the limited lumi- 
nosity ranges of the aforementioned samples. For a given [O HI] lumino sity 
in the high-luminosity, radio-quiet AGN sample of Maiolino et al. 12007), for 
instance, these relations all overestimate the detected X-ray luminosities by 
~l-2 dex. Importantly, applying any of these other correlations would re- 
sult in a significantly larger intrinsic X-ray luminosity estimate, thus favoring 



A08, for instance, have recently shown that several optically- 
identified Compton-thick AGN also follow the AGN correla- 
tion shown in Fig. [2j when one converts uncorrected- [OIII] 
luminosity into X-ray luminosity. We find that the predicted 
rest-frame 2-10 keV luminosities lie a factor of «5-6 above 
the current X-ray limits, supporting the presence of signifi- 
cant X-ray obscuration. Since the optical extinction in the 
xFLS ULIRGs should be much larger than for the type 1 AGN 
used to derive the correlation, these converted X-ray luminos- 
ity constraints should be considered very conservative. 

An alternative scenario sometimes invoked is that the hot 
dust component is somehow associated with massive, deeply 
embedded star formation. Wolf-Rayet and early-type O 
stars, for instance, are certainly capable of producing the 
hard radiation necessary to power the Z^.s^m continua in the- 
ory. Strong continuum-dominated sources like the weak- 
PAH xFLS ULIRGs, however, lie 1-2 dex above the the 
^5.8/jm-ipAH(7.7^m) correlation for star-forming galaxies. This 
effectively rules out star-formation for all but the most con- 
trive d scenarios, but i s fully consistent with known AGN (e.g., 
S07: lDesai et al.ll2007l) . Further indirect evidence against star 
formation comes from high-spatial resolution mid-IR imaging 
of the most powerful local infrared luminous AGN, including 
a few from Table [T] Such observations constrain a large frac- 
tion of the mid-IR continuum light within a region <100 pc in 
size (<0"3) around the nucleus, with mid-IR surface bright- 
nesses in excess of «10 14 L Q kpc" 2 (e.g.. lSoiferetail 2000. 
2003). While the most extreme 'super' star-clusters known 
locally can s till generate compar able mid-IR surface bright- 
nesses (e.g., iGoriian et al.l 1200 ll) . the largest sizes of such 
star-clusters are only 1-1 Opc, which are much smaller than 
the current spatial limits on local ULIRG nuclei. Coupling 
this 1-2 dex size deficit with the facts that the xFLS ULIRGs 
have (1) mid-IR luminosities 1-2 dex higher than the above 
local AGN and (2) relatively weak PAH features, and it seems 
highly unlikely that super star-clusters (even exotic ones at 
high redshift) are capable of producing the required mid-IR 
luminosities of order 10 12 -10 13 L©. As such, AGN appear 
to be the most plausible option for powering the bulk of the 
mid-IR emission in ULIRGs. 

In summary, the strong 5.8 /im continuum and [OIII] lu- 
minosities (once corrected for extinction), coupled with the 
weak X-ray emission even in stacked X-ray images, argues 
for at least mild Compton-thick obscuration in a large frac- 
tion of, if not all, xFLS ULIRGs. 

4.3. Comparison to Other Selection Techniques 

There are several different techniques to find mid-IR AGN, 
and the common goal of these methods is to identify heavily- 
obscured, and potentially Compton-thick AGN that might be 
missed by X-ray, UV, or radio selection. These techniques 
generally fall into a few overlapping categories: (1) those 
like Y07 that look for strong hot dust components relative 
to stellar continua ("mid-IR excess"), (2) those that look for 
hot dust components directly in the IRAC 3.6-8.0/im bands 
alone ("mid-IR spectral slope"), and (3) those that look for 
strong radio emission above that exp ected from sta r form ation 
("radio excess"). A recent study by iDonley et al.l (120081) has 
investigated the reliabilities of these methods. Given the de- 
gree of high quality follow-up data on our sample, particularly 
spectroscopic redshifts, and the likelihood that xFLS ULIRGs 

even more X-ray absorption. 
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FIG. 3. — Rest-frame Ls g^m versus redshift for a select few samples of candidate Compton-thick AGN including: mid-IR excess objects studied here 
(red triangles); X-ray (open green squares) and IR-selected (green circles) objects from P06; radio-excess objects from MS06 (blue triangles); mid-IR excess 
objects from F08/F09 (grey circles for GOODS selection, magenta circles for COSMOS selection). Triangles denote samples with robust redshifts (>90% 
from spectroscopy), while circles denote less-robust samples with redshifts predominantly from photometry (~10% spectroscopy for P06; ~5% and ~40% 
spectroscopy for F08 GOODS and F09 COSMOS, respectively). Aside from the five X-ray detected Compton-thick candidates from P06, we only plot the subset 
of obscured AGN candidates which are individually X-ray undetected for each sample, under the assumption that an X-ray detection of any kind potentially 
implies something less than Compton-thick obscuration. Dashed lines denote expected L2-10 kcV luminosities adopting the rest-frame 5.8 /^m-to-2-10 keV 
correlation shown in Fig. [2] Rest-frame Ls.s/jm values for P06 and MS06 were estimated from interpolation of best-fit power-law models to the published fluxes, 
incorporating an additional boost of fs 1.2 to account for apparent spectral curvature due to extinction; this average boost was determined empirically from the 
xFLS ULIRG sample, where rest-frame 5.8 fim continuum fluxes have been determined directly from spectral decomposition (S07; S08). This estimation method 
is only appropriate for sources with mid-IR spectra well-represented by power-law models, and could deviate somewhat for objects with strong absorption or 
emission features. It should also be stressed that while the Compton-thick AGN nature for the bulk of xFLS ULIRGs appears robust, the same level of confidence 
cannot be assumed for the other samples, for which the relative fractions of sources dominated by star-formation in the IR or obscured but not Compton-thick are 
generally not well established. 



host Compton-thick AGN, a lower limit to the number of po- 
tential Compton-thick AGN discovered by other methods can 
be assessed by investigating the relative frequency with which 
xFLS-like ULIRGs are selected via these other methods. Such 
comparisons also allow us to place the xFLS ULIRGs within 
the context of these other surveys, partic ularly in regard to 
their relative space densities presented in §4A\ Table [3] pro- 
vides an overview of our findings, which we discuss in detail 
below. We automatically exclude xFLS ULIRGs with uncon- 
strained IRAC colors from affected comparisons. Also, since 
sources with direct X-ray detections are likely to be relatively 
unobscured, we only make comparisons between X-ray unde- 
tected subsets of the samples. 

We begin comparisons with the clos ely related studies 
of Martinez-Sansigr e et alJ (120061 120071 collectively here- 
after MS06) and iPolletta et alJ (l2006T hereafter P06), both 
of whom also focus on similarly powerful obscured AGNs 
(see Fig. [3]). The "radio-excess" technique of MS06 selects 
only high-redshift ULIRGs (also within the xFLS footprint) 
which lie above the ra dio-to-far-infra red correlation for star- 
forming galaxies (e.g., Condon 1992). This technique targets 
the radio-intermediate to radio-loud portion of the obscured 



QSOs population, from which the remaining radio-quiet sub- 
set can, in theory, be loosely extrapolated. Within a given 
redshift bin, the average 5.8 [im luminosity of the MS06 sam- 
ple is roughly a factor of two lower, indicating that MS06 
samples less powerful AGN. Importantly, nearly all of the 
MS06 sample has also been o bserved with the Spitzer IRS 
(Martinez-Sansigre et al. 2008), and display properties very 
similar to the xFLS ULIRGs, such as strong silicate absorp- 
tion (T9.7 ~ 1-2) and spectra dominated by strong hot dust 
continua. Under their criteria, five of the xFLS ULIRGs 
sources should nominally have been selected, although three 
of these appear to have been rejected from that sample due to 
source-blending issues in the IRAC 3.6/im band. The fraction 
of xFLS ULIRGs selected by this method thus lies somewhere 
between 5-10%. While this seems low, extrapolating from 
the fraction of radio-intermediate AGN to the full population 
yields a more substantial overlap. 

P06, on the other hand, use a number of criteria (both "mid- 
IR excess" and "mid-IR spectral slope") to remove sources 
with non-power-law SEDs from their IR-selected sample (see 
Table O, selecting sources with a wide range of SED-types, 
luminosities, and obscuration levels. They find that «35% 
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TABLE 3 

Selection Technique Comparison 



Study Selection Criteria xFLS Recovery 

Yan et al. (2007) S^m > 0.9mjy, a(24,8.0) > 0.5, a(24,R) > 1.0 — 

Lacy et al. (2004) 5 2 4 Mm > 4.4 mjy or S$ o M m > Imjy, r(5. 8,3.6) > -0.1, r(8.0,4.5) > -0.2, r(8.0,4.5) < 0.8[r(5.8, 3.6)] + 0.5 Ri65%* 

Stern et al. (2005) ' [R] < 21.5, S 2 4 M m > l.OmJy, ([5.8] - [8.0]) > -0.07, rj30%* 

([3.6]-[4.5]) > 0.2([5.8]-[8.0])-0.16, ([3.6] - [4.5]) > 2.5([5.8] - [8.0])-2.3 

Cool et al. (2006) [/] < 22, S 24Mnl > l.OmJy, ([3.6]-[4.5]) > -0.1, ([5.8] - [8.0]) > -0.05 ps100%* 

Donley et al (2007) S 2 4 Mm > 80 ^Jy, aiRAC < -0.5, where /„ cx v a , P > 0.1 ^70% 

^^^^^^^T (2006) monotonic flux increase in > 3 mid-IR bands, ~70% 

0-2— TAfim < -1, where /„ cx v a , x„ < 13.2(-a 2 _ 24fl „, - 1) < 20, 
r(3.6,g') > 1.18, r(3.6,r') > 1.11, r(3.6/) > 1.0, SED fitting 

Martmez-Sansigre et al. (2006) S 24(im > 0.3rrjy, S 3 6 „ m < 45/xJy, 0.35mJy < Si 5G Hz < 2mJy «5-10% 

Fiore et al. (Mm S 24(im > 0.55mjy (COSMOS) or4()^Jy (GOODS), r(24,R) > 3.0, [R]-[K] > 4.5 ^40% 

Dev et al. (2008) S 24(jra > 0.3mjy, r(24,R) > 3.0 «50% 



Note. — Col. (J) Selection technique. Col (2) Flux/magnitude cutoffs and color selection criteria. Here a(Ai , A2) = log( ^/f^\\^ ) an d , -^2) = log( ^|^j ). Brackets "[A]" denote 
AB magnitudes. A itself denotes the Spitzer IRAC and MIPS band (given as central wavelength in units of /im) or the standard optical/near-IR band. Col. (3) Fraction of sources selected 
by the xFLS ULIRGs criteria which would also be selected by another technique. The fractions denoted with a '*' have been assessed without the stated Sz4^m selection criteria, since 
these would typically reject the bulk of xFLS ULIRGs outright. 



of their sample are X-ray detected, with only two X-ray- 
detected sources showing evidence for Nu > 10 24 cm" 2 ab- 
sorption. The absorption for the X-ray undetected subset of 
their sample should typically be comparable to or larger than 
the X-ray detected sources, but this absorption distribution has 
yet to be firmly established. The 5.8 /tm luminosity and red- 
shift range of the xFLS sources matches the upper end of the 
P06 sample, and -70% of xFLS ULIRGs should be selected 
by the P06 criteria. The large overlap suggests that a sizable 
number of Compton-thick AGN likely exist within the P06 
sample. How this extends to the lower luminosity sources, 
however, and particularly to the unusually large number that 
lie near the flux limit of their survey, is unclear. 

Applying the mos t commonly u sed IRAC AGN color - 
selection cuts from lLacv et al.1 (120041) and IStern et al.l d2005l) . 
we find that w70% and «30% of the xFLS ULIRGs are se- 
lected, respectively. The low percentage for the latter stems 
from the fact that the majorit y of xFL S ULIR Gs he imme- 
diately to t he right of t he IStern et al.l (12005b region. Ex- 
tending th e IStern et al.l (120051) region for z > 1 objects as 
ICool et all (120061) have done, for example, recovers an im- 
pressive «100% of our sample. A major limitation in all 
of the above studies, however, is the adoption of either shal- 
low optical and/or mid-IR flux cutoffs, which effectively re- 
move all xFLS U LIRGs reg a rdless of color-selection crite- 
ria. In the case of lLacv et al.l ( 120041) . the factor of «5 higher 
24/im cutoff, coupled with the expected strongly declining 
space density evolution for ULIRGs, is likely to severely limit 
the number of comparably luminous but mor e ne arby xFLS- 
like U LIRGs detected. For lStern et all (T2005) and lCool etaTI 
(2006), the bright optical selection will effectively remove a 
large fraction of highly-obscured ULIRGs like those studied 
here. While extending the above techniques to deeper sam- 
ples would help select ULIRGs like those in the xFLS, do- 
ing so could be counterproductive since it coul d also intro- 
duce many new star-forming contaminants (see Donley et al. 
2008). Adding an additional color selection to remove z < 1 .2 
galaxies, such as an R - K criteria, could effectively address 
this apparent limitation. 

The other main "mid-IR spectral slope" technique used 
extensively is the power- law galaxy (PLG) crit eria of 
Alonso-Herrero et all (120061) and IDonley et all (120071) . which 
Donley et al.l d2008) contend recovers the majority of high- 



quality AGN candidates. All of the xFLS ULIRGs would 
qualify as PLGs based on their underlying mid-IR AGN con- 
tinuum slopes as found from spectral decomposition. How- 
ev er, using only the o bserved IRAC data alone as stipulated 
in IDonley et al.l (I20071) . only r;70% of the xFLS ULIRGs 
are recovered as PL Gs. While not in direct contrast with 
IDonley et al.l £2008), our result implies that there exists a 
rather substantial population of legitimate mid-IR excess 
AGN which are currently missed by the PLG technique. 
These non-PLG objects typically either had poor IRAC con- 
straints (and subsequently large uncertainties in their spectral 
slopes) or were dominated by stellar continua well into the 
IRAC bands but had strong hot-dust components at longer 
wavelengths (see S07). Intriguingly, several of the PLG- 
selected xFLS ULIRGs also had relatively strong host stellar 
continuum components, demonstrating that the PLG method 
does select some composite AGN. 

Among th e "mid-IR excess" techni ques, the [Fiore et al.l 
d2008ll2009l hereafter F08 and F09) and bev et al.l d2008l) se- 
lection criteria are the most closely related to Y07, although 
both employ much lower 24/im flux cutoffs, have no 24/im- 
to-8/im flux selection, and require a factor of «3 larger 24/im 
to R-band excess. F08/F09 additionally implement an R-K 
cut to remove low-redshift contamination from star-forming 
galaxies. The lower flux cutoff there probes deeper into the 
overall 24/im source population to select both lower luminos- 
ity and higher redshift objects (see Fig. 0, while the larger 
mid-IR excess picks only the reddest ULIRGs. The lack of 
24/im-to-8/im selection, however, should lead to the inclu- 
sion of many sources with either less obscuration and/or hot- 
ter dust continua. We find that only «40% and «50% of 
xFL S ULIRGs lie at the extremes selected by the F08/F09 
and iDey et all (120081) techniques, respectively. There is a 
mild tendency for these techniques to predominantly select 
higher redshift sources as expected, although several high-z 
xFLS ULIRGs are still notably excluded. If all of the 5.8/im 
luminosity in F08/F09-selected objects stems from obscured 
AGN activity, then the X-ray stacking detection performed 
by these authors equates to an average mid-IR/X-ray decre- 
ment of w 800, which is an order of magnitude below our 
xFLS X-ray luminosity stacking upper limits, and hence fully 
consistent with o ur results. Anothe r mid-IR excess selection 
method is that of iDaddi etail (120071 hereafter D07), which 
selects a large number of high-redshift, moderate-luminosity, 
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FIG. 4. — Co-moving space density of candidate Compton-thick 
AGN in three approximate rest-frame, intrinsic X-ray luminosity ranges: 
i<2-l0feV = lO^-lO 44 erg s _1 (squares), L2-l0feV = 10 44 — 10 45 erg s (cir- 
cles), and L2-iokeV = 10 45 -10 46 erg s (stars). Only a subset of constraints 
from Table \4\ are shown. The color-coding of the samples is identical 
to Fig. [J] with the additions of Compton-thick candidate AGN within the 
GOODS region from D07 (purple) and A08 (orange). Similar to Fig. [J] 
samples with mostly robust spectroscopic redshifts are shown as filled sym- 
bols, while less-robust photometrically-constrained samples are shown with 
open symbols. For comparison, we show the predicted space densities of 
Compton-thick AGNs from the G07 CXRB model for these three X-ray lu- 
minosity bins (dashed lines). The number of candidate Compton-thick xFLS 
ULIRGs are roughly equal to the predicted values to within errors in the 
i-2-lOfev = 10 44 -10 45 erg s -1 bin, while P06 and F09 lie a factor of ~3^1 
higher. The fact that these sample selection techniques do not fully overlap 
implies that Compton-thick QSOs are likely ubiquitous at high redshift. 



obscured AGN candidates comparable to the GOODS sources 
in F08. This selection technique is too complex to employ 
here, since it relies on assumption-dependent excesses be- 
tween the dust-corrected UV and IR star formation rates in 
each object. Nonetheless we note that the D07 sample has an 
average 5.8/im excess luminosity of ~3 x 10 44 erg s -1 , which 
is comparable to those of our z~ 1 objects, and an X-ray stack- 
ing detection an order of magnitude below our z ~ 1 X-ray 
luminosity upper limits, which again is fully compatible with 
our results. 

While it is rare to find full consistency between Y07 and 
other selection methods, there appears to be enough overlap 
(if we neglect shallow flux criteria) to demonstrate that many 
of these methods select a substantial subset of powerful, po- 
tentially Compton-thick ULIRGs like the ones characterized 
here. Aside from MS06, however, a fundamental limitation of 
these other studies is their reliance on broad-band photometry 
and, in several cases, photometric redshifts. This critical lack 
of precise redshifts and spectrally-deconvolved SEDs makes it 
impossible to distinguish unambiguously between mid-IR ex- 
cesses due to hot dust from AGN and those from strong PAH 
features due to vigorous star formation. While this may be 
a relatively small concern at high mid-IR luminosities where 
AGN are thought to dominate (e.g., Ls.^m^lO 45 erg s -1 ), it is 
likely to cause severe problems at lower mid-IR luminosities 
where star formation is expected to be u biquitous (e.g., S07; 
lLacv et al.ll2007l: iMartmez-Sansigre et al.ll2008t) . 

4.4. 77ie Space Density of Heavily-Obscured Accretion 



We now investigate the space densities of the candidate 
Compton-thick AGN samples discussed in ^4.3l in several red- 
shift and predicted intrinsic X-ray luminosity bins. Approx- 
imate values of intrinsic rest-frame L2-iokeV were predicted 
from rest-frame L^.s^m assuming the conversion above. 16 A 
key concern is that although these samples all span a wide 
range of luminosities and redshifts (see Fig.[3]l, their selection 
is often potentially biased or incomplete within a particular 
bin due to limited solid angle coverage, restrictive selection 
criteria, or survey flux limits. Table [Uprovides space density 
estimates for the full range of parameter space probed, while 
Fig. E] highlights the best constraints for each of the samples. 
For bins which are obviously incomplete, we assume the de- 
rived values are lower limits. We again strongly caution over 
interpretation of sources in the comparison samples which of- 
ten lack both complete optical/near-IR spectroscopic identifi- 
cation and well-sampled SEDs with IRS spectroscopy to de- 
couple possible emission mechanisms. Contamination from 
star formatio n, pa rticularly at lower mid-IR luminosities (as 
discussed in ^4. 3b . could be substantial. Since this contam- 
ination is not yet well-constrained, we liberally place these 
comparison samples at their highest space densities using all 
candidate objects. 

For comparison, we also show in Fig. @] the CXRB syn- 
thesis predictions from G07 for three different luminosity 
ranges of Compton-thick AGN. These curves should be re- 
garded as approximate, since considerable uncertainty re- 
mains in the redshift, luminosity, and column density dis- 
tributions of the sources that comprise the full CXRB. In 
particular, high-luminosity AGN such as these considered 
here are predicted to contribute relatively little power to 
the overall CXRB, and conversely cannot be strongly con- 
strained by such models. We note that G07 adopts a model 
for intrinsic rest-frame AGN with L2-10 keV^lO 44 ergs s" 1 
consisting of equal parts unobscured (Nh < 10 22 cm -2 ), ob- 
scured Compton-thin (Nh ~ 10 22 -10 24 cm" 2 ), and Compton- 
thick AGN. They consider all AGN with N H > 10 22 cm -2 
(both Compton-thin and Compton-thick) as obscured and 
Nh < 10 22 cm -2 as unobscured, such that their model implies 
an expected obscured-to-unobscured ratio (hereafter simply 
"obscured fraction") of 2:1 among luminous AGN with the 
G07 model. We stress that the space densities of unobs cured 
QSOs (e.g.. lHasinger et alj 120051: iHqpkins et all 120071) and 
Compton-thin QSOs (e.g., lUeda et aljf2003t La Fran ca et alj 
l2005t ISilverman et alj 120081) are well-constrained by obser- 
vations, while the number of Compton-thick AGN have only 
been inferred from modeling. We investigate here how many 
Compton-thick AGN might be expected from several mid-IR 
selected surveys. Because these three classes of AGN are 
equal in the G07 model, the lowest two dashed Compton- 
thick AGN curves also represent the expected space densities 
of unobscured and obscured Compton-thin QSOs; as such, the 
lines serve as useful visual benchmarks to assess the potential 
obscured fractions among various samples. 

Considering first the AGN between predicted 

L 2-io kev~ 1()45 - 1046 er S s s ~'' we find that the xFLS ULIRG 
sample yields a space density of $ ~(1.44^4g)x 10" 7 Mpc" 3 
at z = 2-3, where the xFLS selection technique was designed 



16 The choice of correlation slope and intercept do not have a strong effect 
on the estimated luminosities (and there space densities), and lead to varia- 
tions of perhaps f»2 in either direction. Additionally, correcting rest-frame 
i-5.8^m for extinction could raise it a factor of Ri2 on average, implying a 
similar increase to intrinsic rest-frame Z-2-lOkeV- 
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TABLE 4 

Compton-Thick Candidate Space Densities 



Sample 


Area 


z 


i2-l0teV = 10 44 -10 43 erg s 


La-\0keV = 10 43 -10 4b erg s 


Yan et al. 2007 (xFLS) 


3.26* 


0.5- 


-1.0 


<l.44E-07(0) 


<1.44E-07 (0) 






1.0-2.0 


>l.32+°^E-07 (18) 


<1.64E-07 (0) 






2.0-3.0 


>1.85+jj;^E-07(9) 


1.44+°-«|E-07 (7) 






3.0^1.0 




>2.16+ 4 |?E-08(1) 


Yan et al. 2007 (xFLS) 


3.26* 


0.5- 


-1.0 


<1.44E-07 (0) 


<1.44E-07 (0) 


w/o strong PAH sources 




1.0-2.0 


>9.53 + |'^E-08 (13) 


<1.64E-07 (0) 






2.0- 


-3.0 


>1.44° 6 8 6 E-07(7) 


1.03+jj ^E-07 ( 5 ) 






3.0^1.0 


— 


>2.16+^E-08 (1) 


Martmez-Sansigre et al. 2006 (xFLS) 


3.70 


0.5- 


1.0 


<5.08E-07 ( 0) 


<5.08E-07 (0) 




1.0- 


-2.0 


>3.10 + :J : ^E-07 (12) 


<1.45E-07 (0) 






2.0- 


-3.0 


>1.36^'^E-07 ( 6) 


2.26+? S5E-08 (1) 






3.0^1.0 


— 


>2.38+-; : ^E-08 (1) 


Polletta et al. 2006 (SWIRE) 


0.60 


0.5- 
1.0- 


-1.0 
-2.0 


2.24+' o,E-06 ( 4) 
>4.15+° : ™E-06(26) 


<3.13E-06 (0) 
<8.93E-07 (0) 






2.0- 


-3.0 


>2.93+°;™E-06 (21) 


5.58+| gE-07 (4) 






3.0^1.0 


>2.93 , sg E-U/ ( I) 




Fiore et al. 2009 (COSMOS) 


0.90 


0.5- 


-1.0 


1.12 + ''»E-06(3) 


<2.09E-06 (0) 






1.0- 


-2.0 


>3.62+° jgE-06 (34) 


2.13+^E-07 (2) 






2.0-3.0 


>1.49^E-06(16) 


6.50+| : ^E-07 (7) 






3.0-4.0 




>2.93+f : ^E-07 (3) 


Fiore et al. 2008 (GOODS-S) 


0.04 


0.5- 


-1.0 


<4.73E-05 ( 0) 


<4.73E-05 (0) 






1.0- 


-2.0 


2.65+'!"E-05 (11) 


<1.35E-05 (0) 






2.0- 


-3.0 


101 -0 15 E_04 ( 48 > 


<1.18E-05 (0) 






3.0-1.0 


1.99+[}^E-05 (9) 


<1.24E-05 (0) 


Alexander et al. 2008 (GOODS-N) 


0.04 


2.0- 


-2.5 


1.50+^E-05(4) 





NOTE. — Col. (1) Sample. Col. (2) Area of sample, in units of deg~ . '*' — Note for our sample, the original area (3.7 deg~ ) has 
been reduced by the fraction of original sources selected by the Y07 criteria (59) but not followed-up with Spitzer IRS; see Y07 for further 
details. Col. (3) Redshift range. Cols. (4-5) Space densities derived from number of sources (in parenthesis) in the rest-frame 2-10 keV 
luminosity ran ges of 10 44 -10 4 5 erg s _1 and 10 45 -10 46 erg s~ , respectively, as shown in Fig.f3] Errors quoted are based on counting statistics 
only assuming Gehrels i 1986). and are quoted at la confidence when a source exists within a bin and 3a for upper limits. Bins shown as 
lower limits suffer from incompleteness due to flux limits of the various surveys, as can be seen in Fig. [3] 



to be most sensitive. This value is already «90% of that 
predicted by the G07 model. However, as determined in 
^4.21 under the "conservative" A^h > 10 24 cm" 2 selection 
criteria, assuming the IRS-derived 5.8 fj,m AGN continuum 
luminosities provide a good proxy for the expected X-ray 
luminosity, our subset sampling statistics only allow us 
to cautiously assume that >25% of the 0.5 < z < 2.0 
and >80% of 2.0 < z < 4.0 xFLS ULIRGs here are 
Compton-thick at 90% confidence. Assuming equal parts 
unobscured and Compton-thin AGN, the above constraint 
on the Compton-thick space density becomes >70% that 
of G07, or equivalently an obscured fraction of > 1 .7: 1 

among powerful L2-10 kev^^ 45- ^ 46 er § s s_1 QSOs. Our 
other high-luminosity redshift bin is a factor of a few 
lower, but remains consistent with the z = 2-3 prediction 
due to potential incompleteness. In the lower range of 

^2-10 keV^lO 44- ^ 45 er S s s_1 > however, the xFLS detects 
only ~10% of the objects predicted by G07, illustrating 
that the 24/im flux limit is far too shallow to detect the 
majority of these potential faint, high-z ULIRGs; note that 

the ^2-10 keV^lO 44- !^ 45 er 8 s s ~' vaiues from the xFLS and 
MS06 samples have not been plotted in Fig. [4] to reduce 
visual clutter, but are listed in Table|4]for completeness. 

As shown in Table [4] excluding strong-PAH ULIRGs 
where the Ls^m continuum measurements are not as well- 
constrained does not strongly affect our results. For instance, 
even if we pessimistically neglect the two strong-PAH sources 
in the z = 2-3 bin, the obscured fraction is still > 1.5:1. We 
caution that these estimates reflect counting and sampling 



statistics for the xFLS sample only, and do not account for 
systematic errors such as the slope or dispersion in the Ls.s^m- 
^2-io keV conversion, the unknown mid-IR extinction correc- 
tions, or sample selection completeness for Compton-thick 
AGN. As we have argued in previous sections, these ef- 
fects largely skew toward substantially underestimating the 
true number of Compton-thick AGN and the Ls.s^m values 
of such AGN, ensuring that our derived obscured fraction 
limits should be relatively robust. As such, our results im- 
ply that the simple mid-IR selection criterion laid out in Y07 
can be highly efficient for finding large numbers of powerful 
Compton-thick QSOs candidates. To pick up more modest 
(and typical) Compton-thick QSOs, however, we must move 
to deeper mid-IR surveys. 

Looking at the estimates from the P06 and F09 techniques, 
both appear to find exceptionally high numbers of obscured 
QSO candidates in several redshift and predicted luminosity 
ranges. At ZJJ-io kev~10 45 -10 46 er g s s > these estimates lie 
a factor of «3 above the G07 curve, implying a huge pop- 
ulation of heavily-obscured AGN and an obscured fraction 

of perhaps >4:1. Even in the ^2-10 kev~^ 44- ^ 45 er 8 s s ~' 
regime, where these surveys become incomplete, these stud- 
ies still find space density lower limits at roughly comparable 
levels to the G07 predictions. Intriguingly, these constraints 
for both luminosity ranges are consistently high across several 
adjacent redshift bins, indirectly arguing against strong con- 
tamination by wide PAH features, which we might expect to 
dominate primarily around z ~ 2 when rest-frame PAH lines 
pass through the 24/im bandpass. 

Among the deeper pencil-beam surveys, we show the 
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results of A08 and F08 in the GOODS regions for 

L 2-io kev~ 1044 - 1045 er S s s_1 AGN - The A08 sample consists 
of four relatively secure Compton-thick AGN candidates in 
the GOODS-N region at z = 2-2.5, all of which have robust 
optical/mid-IR spectroscopy and deep X-ray constraints. The 
space density of such sources is a factor of sa2 higher than 
the G07 model curve and implies an obscured fraction of 
perhaps «3:1 overall. The F08 sample, on the other hand, 
employs the same method as F09 but to substantially fainter 
24/im fluxes. As such, they find Compton-thick AGN candi- 
date space densities as high as a factor of ^4-20 above the 
G07 predictions, implying an fraction well in excess of «5:1. 
In contrast to the consistency across redshift bins for the 

ZiLio keV^lO 45- !^ 46 er S s S_I sources selected by this method, 
however, the z = 2-3 bin here is a factor of 5 higher than adja- 
cent bins. While this could be a legitimate evolutionary trend, 
such a high outlier arouses suspicion and may perhaps best be 
explai ned by contaminati on from star-forming PAH features 
(e.g.. iMurphy et afl 12009). 17 For comparison we also show 
the constraints of D07, who probe slightly lower AGN lumi- 
nosities than F08, in the range L^ lQkeV mlQ 43 -10 44 ergs s" 1 
over the redshift range z= 1.4-2. 5. Their space density con- 
straint is similarly high, with no data to compare against at 
higher or lower redshifts. Like F08, this dataset also pre- 
sumably suffers from uncertainties both in AGN luminosity 
determinations and in the fraction of the sample which con- 
tribute to the AGN signal. The disparity of the F08 and D07 
results in the z ~ 2 regime highlights the need for deep mid- 
IR spectroscopy to dispel the current ambiguity among the 
fainter candidate Compton-thick AGN population. 

In general, it appears that the relatively conservative xFLS 
criteria selects a substantial number of strong Compton-thick 
AGN candidates which are consistent with current predic- 
tions, but should be considered by no means exhaustive. Other 
techniques such as those adopted by F09 or P06 could present 
immensely effective methods for identifying obscured AGN 
population, at least for high-luminosity samples. Clearly sub- 
stantial caution must be exercised until individual AGN and 
star formation contributions can be effectively pinned down 
in these samples. Nonetheless, even if only a minority of the 
candidates are eventually confirmed, the sheer number of lu- 
minous Compton-thick AGN candidates found by all of the 
various selection techniques, which only partially overlap, 
argues in favor of much higher obscured fractions at these 
extreme AGN luminosities than have been determined from 
X-ray-detected studies alone. Such fractions could be well in 
excess of 4: 1 if these various techniques indeed select a high 
percentage of secure Compton-thick AGN. Where this leaves 
recent claims about th e luminosity-dependence of obscuration 
(e.g.,|Hasinger 2008) is unclear. 

5. CONCLUSION 

The Chandra X-ray upper limits we report here indicate 
that the substantial obscured AGN accretion found among 
mid-IR selected, high-redshift ULIRGs in the xFLS is likely 
to be Compton-thick in nature. The number density of 
these candidate Compton-thick AGN already nearly equals 
that of unobscured QSOs and obscured Compton-thin AGN 
(e.g., such as the co mpa rably numerous A GN studied by 
Weedmanetal. 2006] and [BrandetaLl [2008), suggesting an 



obscured fraction of at least w 1.7:1. We regard this con- 
straint as a robust lower limit, since the assumptions made 
to achieve this result tend to substantially underestimate 1) 
the true number of Compton-thick AGN since our mid-IR se- 
lection only partially overlaps with candidate Compton-thick 
AGN selected by other techniques, and 2) the Ls.s^m values of 
the AGN we do select. This fraction could easily exceed 4: 1 
among QSO-luminosity AGN if even a fraction of the candi- 
dates from various selection methods are confirmed by mid- 
IR spectroscopy. 

We stress that the uniqueness of our sample, with its rel- 
atively simple selection criteria, complete Spitzer-IRS spec- 
troscopy, and UV-to-radio SEDs, provides several indepen- 
dent lines of evidences for powerful AGN activity, and should 
be regarded as far more reliable than similar photometric- 
based samples. This is a major concern particularly for fainter 
mid-IR samples (e.g., D07; F08) which attempt to tease out 
small AGN contributions from large samples of objects with 
dominant starburst components. Unfortunately, the neces- 
sary mid-IR spectroscopic identification and deconvolution 
required to confirm these controversial faint sample results 
must await the launch of the James Webb Space Telescope. 

Such large obscured AGN fractions at high luminosities 
have a few important implications. First, this result adds 
uncertainty to recent optical and X-ray constraints that call 
for some luminosity-dependence in obscuration and/or the 
opening angle of the p utative AGN torus (so-called "reced- 
ing torus" models. e.g.. | Lawrencell9'9ll Hasinger et al.ll200"5l 
Maiol ino et"aTI l2007t lHasingeri 120081) . Whether this appar- 
ent discrepancy in number density ultimately stems from a 
short (but obviously ubiquitous) evolutionary stage of ob- 
scured growth or a more fundamental, long-standing prop- 
erty of the immediate AGN environment remains to be seen. 
If similar obscured fractions can be confirmed over a larger 
range of redshifts and luminosities, it would certainly place 
significantly stronger constraints on the form and nature of 
SMBH evolution, particularly in relation to satisfying both 
the local bulge luminosity and CXRB constraints. Sec- 
ond, the implied number of such powerful, dust-laden AGN 
may provide a plentiful driver for AGN feedback scenarios 
(e.g.. ISilk & Reesl [19981: lHaehnelt et all [T998I: iFabianl [19991: 
iMurrav et al.l 120051: iFabian et alj |2008) and may ultimately 
lead to considerable enrichment of intergalactic material (e.g . , 
iDave et al.l200ltlAdelberger et alJl2003llMenard et alT2009l) . 
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17 This many heavily-obscured AGN could also be difficult to rec- |2005tlShankar et alj2007l) 
oncile with limits imposed by the local black hole ma s s density (e.g., 
Yu&Tremaine 2002; Marconi et al. 2004; Merloni 2004; La Franca et al. 
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